Journal of Low Temperature Physics manuscript No. 

(will be inserted by the editor) 



L. Gironi 1 



Pulse Shape Analysis with scintillating 
bolometers 



01.08.2011 



Keywords Bolometers, Scintillators, Pulse shape discrimination 

Abstract Among the detectors used for rare event searches, such as neutrinoless 
Double Beta Decay (OvDBD) and Dark Matter experiments, bolometers are very 
promising because of their favorable properties (excellent energy resolution, high 
detector efficiency, a wide choice of different materials used as absorber, ...). How- 
ever, up to now, the actual interesting possibility to identify the interacting parti- 
cle, and thus to greatly reduce the background, can be fulfilled only with a double 
read-out (i.e. the simultaneous and independent read out of heat and scintillation 
light or heat and ionization). This double read-out could greatly complicate the 
assembly of a huge, multi-detector array, such as CUORE and EURECA. The 
possibility to recognize the interacting particle through the shape of the thermal 
pulse is then clearly a very interesting opportunity. 

While detailed analyses of the signal time development in purely thermal de- 
tectors have not produced so far interesting results, similar analyses on macro- 
bolometers (^10-500 g) built with scintillating crystals showed that it is possible 
to distinguish between an electron or y-ray and an a particle interaction (i.e. the 
main source of background for OvDBD experiments based on the bolometric tech- 
nique). Results on pulse shape analysis of a CaMo04 crystal operated as bolome- 
ter is reported as an example. An explanation of this behavior, based on the energy 
partition in the heat and scintillation channels, is also presented. 
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1 Introduction 

In neutrinoless Double Beta Decay (OvDBD) 1 , as in all the rare event studies, 
spurious events are a limiting factor to the reachable sensitivity of the experiment. 
Unfortunately natural radioactive background is often present in the detector itself 
or in the materials surrounding it, no matter how much one can try to reduce it 
with shieldings. In order to reduce this residual unavoidable background, all the 
possible approaches require both a good energy resolution and the capability to 
identify the nature of the projectile that interacted with the detector. 

Bolometers 2 are based on the detection of phonons produced after an energy 
release by an interacting particle and can have both an excellent energy resolu- 
tions and extremely low energy thresholds with respect to conventional detectors. 
They can be fabricated from a wide variety of materials, provided they have a low 
enough heat capacity at low temperatures. The latter is a priceless feature for ex- 
periments that aim at detectors containing particular atomic or nuclear species to 
optimize the detection efficiency. If other excitations (such as ionization charge 
carriers or scintillation photons) are collected in addition to phonons, bolometers 
have already shown to be able to discriminate nuclear recoils from electron recoils, 
or a particles from j3 particles and y-rays. Recently has been demonstrated— the 
possibility to obtain similar results by Pulse Shape Analysis (PSA), without the 
requirement of a double readout for phonons and ionization or scintillation light. 



2 Scintillating Bolometers and Pulse Shape Analysis 

The concept of a scintillating bolometer is very simple: a bolometer coupled to a 
light detector 4 . The first must consist of a scintillating absorber thermally linked to 
a phonon sensor while the latter can be any device able to measure the emitted pho- 
tons. Excellent results were obtained using as light detector a second 'dark' thin 
bolometer 5 . In this device the simultaneous and independent read out of the heat 
and the scintillation light permits to discriminate events due to j3/y, a and neu- 
trons thanks to their different scintillation yield. Dark Matter as well as OvDBD 
searches can benefit of this capability of tagging the different particles, and more 
generally this technique can be exploited in any research where background sup- 
pression or identification is important. 

The development of such a hybrid detector, able to discriminate a particles, 
was the main purpose of studies on scintillating bolometers for OvDBD searches. 
Several devices 6 , differing mainly in the scintillating crystal material and size, 
were tested to study their thermal response, light yield and radio purity. In these 
studies was discovered an extremely interesting feature of some of the tested crys- 
tals: the different pulse shape of the thermal signals produced by a's and by j3/y's. 
This feature opens the possibility of realizing abolometric experiment that can dis- 
criminate among different particles, without the need of a light detector coupled 
to each bolometer. In the case of a huge, multi-detector array, such as CUORE 7 
and EURECA^, the benefits of employing this technique would be impressive: 

- more ease during assembly because the single element of the array would be a 
quite simpler device and light collectors are no more needed. 
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- fewer readout channels, with an evident reduction of cost and work and also a 
cryogenic benefit: a significant reduction of the heat load of the readout chan- 
nels. 

- a significant cost reduction, saving money and work that would be necessary 
for the light detectors procurement and optimization. 

Up to now, the discrimination based on PSA has been very effective for par- 
ticles with energy higher than few hundreds of keV. This technique can therefore 
already be applied to OvDBD experiments while can not yet be applied to study 
low-energy events as required for Dark Matter experiments. However, since these 
are the first observations of this feature, there are considerable rooms for improve- 
ment, for example by using phonon sensors and electronics better suited to this 
type of measurements than those used (Ge-NTD thermistors). This could signifi- 
cantly lower the energy threshold of the pulse shape discrimination. 

Finally, the removal of the light guide - used in scintillating bolometer to 
collect more light - allows to face these devices to the assembly materials by 
opening the possibility to study surface contaminations of these materials with 
higher sensitivity than standard technique (e.g. Si surface barrier detectors). The 
study of such contamination is very interesting because they are indicated as the 
main source of background for OvDBD experiments based on the bolometric tech- 
nique^. 



3 Experimental Results 

Several experimental evidences 3 have shown the possibility to discriminate the 
type of interacting particles through different shape of the thermal signal measured 
in different scintillating materials (CaMoO/t, MgMoO/t, ZnMoO/t, ZnSe). As an 
example, I will report here results obtained with a CaMoCU crystal with a mass 
of 158 g (h = 40mm, = 35mm) operated as bolometer in a dilution ( 3 He/ 4 He) 
refrigerator installed in the hall C of the Underground Laboratory of Gran Sasso. 

In figure Q] is reported the thermal rise time T„> (evaluated on the recorded 
raw-pulse as (t9o%-tio%)) as a function of the energy released in the bolometer. 
Two separate bands, ascribed to /3/y and a events can be identified. With differ- 
ent color are reported /3/y and a events having an energy between 2 and 4 MeV 
where are expected to attend the peaks of the most promising OvDBD isotopes. A 
gaussian fit of these two distribution yields a rise-time of (5.788±0.017) ms for 
/3/y's and of (5.649±0.013) ms for a's resulting in a discrimination confidence 
leveQ of Dffi Jf ,=6.5a. 

Similar results were obtained with other crystals (MgMo04, ZnMoO/t, ZnSe) 
showing that this is a common feature of some scintillating crystals. Instead, up to 
now, no evidence of this feature was observed in purely thermal detector. 

The different shape of thermal signal produced by a particles and f}/y particles 
could be explained by the dependence of light yield on the nature of the interacting 
particle. The high ionization density of a particles implies that all the scintillation 

1 The discrimination confidence level D#, Jf is defined as DRj se = R ' se ^/j' Rlsea where Risen i y 

and Rise a are the average rise times produced by the two kinds of particles and dp iy and a a is 
the width of the two distributions. 
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Fig. 1 (Color online) Scatter plot of t,-, jc vs. Heat in CaMoCU crystal. In green (light in the 
upper curve) events in the 2-4 MeV region due to /3/y particles and in red a events. The energy 
calibration is performed on the y peaks. The total live time of this measurement was ~43 h. 
In order to have a high number of a counts in the 2-4 MeV region in these measurements a 
degraded 238 U source— was placed in front of the crystal. 

states along their path are occupied 10 . This saturation effect does not occur or at 
least is much less for /3/y particles. Therefore, in a interactions a larger fraction of 
energy flow in the heat channel with respect to j3/y events. This leads not only to a 
different light and heat yield but also to a different time evolution of both signals. 
The pulse shape of the thermal signal then can be explained by the partition of 
energy in the two channels. 



4 The energy partition and the PSA 

A possible interpretation of the pulse shape discrimination observed in some scin- 
tillating bolometers is based on the following assumption: the processes that lead 
to the thermal pulse formation (i.e. the production of scintillation light and phonons 
and their collection in the sensor) can not be considered instantaneous even in slow 
detectors such as the bolometers. Indeed, at low temperature the averaged scintil- 
lation decay time of molybdates is of the order of hundreds of /is and there- 
fore it is comparable to the typical rise time of the heat signal of our scintillating 
bolometers. 

In the following a possible simplified model to explain the pulse shape dis- 
crimination observed in scintillating bolometers will be presented. Let us assume 
that the interaction of a /3 particle with the crystal produces Nj3 electron/hole pairs. 
A fraction of these ((Op) will recombine through scintillation processes while the 
remaining (l-(Op) will give rise to thermal processes (this holds in the case in 
which only two channels share the total energy). For an a particle the fraction 
that will recombine through scintillation processes ((O a ) will be lower ((O a < (Op) 
because of saturation effects. In fact, most scintillating centers are occupied and a 
larger fraction of energy flow through the heat channel. 
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Since the thermal processes have a very fast time evolution (< 10 s) they can 
be considered instantaneous with respect to thermal pulse in macro-bolometers. 
Therefore, let us analyze more in detail the scintillation channel. 

In figure |2] is shown a schematic representation of the scintillation channel. It 
could be divided in 3 stages. In the first stage (B-B'), which includes from the 
initial energy deposition to the formation of excited luminescent centers, there is 
the emission of 'prompt' phonons with typical time < 10 _8 s^. These phonons 
are then emitted with those of the heat channel and do not give any information 
regarding the type of particle that has released energy. The first stage is followed 
by the emission of scintillation light (B'-A')- This process is very slow if compared 
to the production of phonons. For instance the measured scintillation decay time 
constant for CaMoCU is 380 ±20 /is at 7K 14 . The scintillation process is then 
followed by a further stage of non-radiative de-excitation (A'-A). Phonons emitted 
from this stage are then produced at a later time (i.e. they are produced with an 
exponential behavior equal to the decay time of scintillation) with respect to the 
'prompt' phonons. 




Fig. 2 Schematic configuration coordinate diagram of the scintillation process. The energy E is 
plotted versus the coordinate Q (configurational coordinate in the lattice). 



The fraction Tj of the scintillation channel that gives rise to non-radiative pro- 
cesses is not negligible. For some molybdate s 14 i 15 it could be as high as ^50% of 
the energy that flow in the scintillation channel. Taking into account that the light 
yield of molybdates at low temperature is higher than tens of keV/MeV— , this 
fraction could modify the total thermal pulse shape. In summary the development 
of the total thermal signal caused by the interaction of a j3 particle will be due the 
sum of the heat channel ((Op ) and the non-radiative part of the scintillation channel 
T](l - (Op ) that is partly emitted later. 

In the case of the interaction of an a particle there will be an higher percentage 
of non-radiative processes ((O a > (Op) resulting in a different time evolution of the 
total thermal process. The rise and the decay time for a events measured with the 
CaMo04 crystal are <-^200 /is faster than those of /3/y events. 

Based on this model, the observation of this effect is therefore related to the 
combination of two factors: the scintillation decay time (which should be long 
enough) and the percentage of non-radiative de-excitation of the scintillation chan- 
nel. In addition, the processes of self-absorption should add up to the effect just 
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described. Therefore, it is expected a dependence on the shape and size of the 
crystal. 

Finally, the processes described above relate only to the energy subdivision 
and to the time evolution of the heat and light channels. The heat channel induces 
the temperature rise of the scintillating crystal while the restoration of the equilib- 
rium temperature is mainly due to the thermal link to the heat bath. Further studies 
are necessary in order to include this simplified model in a more general thermal 
model of scintillating bolometer that could be able to describe the measured ther- 
mal pulse shape. 



5 Conclusion 

The possibility to discriminate the nature of the particle interacting in a bolometric 
detector, simply on the base of the shape of the thermal pulse, is proved and opens 
new possibilities for the application of these devices in the field of rare events 
searches. A possible explanation of this behavior based on the energy partition 
in two channels (heat and light) is presented. Further studies are therefore being 
developed to analyze the convolution of the phonons production process and the 
detector thermal response. 
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